Superconducting double spin valve with extraordinary large tunable 

magnetoresistance 
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A superconducting double spin valve device is proposed. Its operation takes advantage of the 
interplay between the spin-filtering effect of ferromagnetic insulators and superconductivity-induced 
out-of-equilibrium transport. Depending on the degree of nonequilibrium, extraordinary large tunnel 
magnetoresistance as large as 10 2 ...10 6 % can be obtained for realistic material parameters, and it 
can be tuned over several orders of magnitude under proper voltage biasing and temperature. The 
relevance of this setup for low-temperature applications is further discussed. 

PACS numbers: 72.25.-b,85.75.-d,74.50.+r 



The possibility to generate highly spin-polarized elec- 
tric currents is one the main topics in the field of spin- 
tronics [l[, both from the fundamental and the tech- 
nological point of view [J]. In this context, a widespread 
route to achieve spin-polarized transport is through tun- 
neling magnetoresistance (TMR) [1| in magnetic tun- 
nel junctions, which consist typically of ferromagnetic 
electrodes separated by a nonmagnetic insulating layer. 
Magnetic tunnel junctions have recently demonstrated 
TMR values as large as 500% at room temperature @, 
and show nowdays an enormous technological impact for 
the realization of a wide range of devices spanning, for in- 
stance, from magnetic random access memories (MRAM) 
to sensors. On the other hand, tunneling through a ferro- 
magnetic insulator (FI) @, 0, Q from nonmagnetic elec- 
trodes is an alternative and effective choice to produce 
highly spin-polarized currents. Here one takes advan- 
tage of the spin-filtering property of FIs due to the dif- 
ferent barrier heights experienced by the two spin species 
(9H10I 11, 12 1 . Spin filtering efficiency approaching 100% 
has been achieved in tunnel junctions comprising Eu- 
based FIs [HI EH, 0, EH, thus showing the potential 
of these materials for spintronic applications. 

Here we propose a double spin valve device 
which combines the spin-filtering effect of FIs with 
superconductivity-induced nonequilibrium. Extraordi- 
nary large TMR (~ 10 2 ...10 6 %) can be achieved, which 
can be tuned over several orders of magnitude under volt- 
age biasing and temperature. Such features make this 
structure attractive for the implementation of magne- 
toresistive nanodevices in superconducting spintronics. 

The system we envision consists of two equal supercon- 
ductors (S) symmetrically connected to a normal-metal 
layer of thickness tw through ferromagnetic insulating 
barriers (FI) [see Fig. 1(a)]. We assume the FI layers to 
be identical (a barrier asymmetry would not change the 
overall physical picture) , d labels the barriers width, and 
a voltage V is applied across the structure. The mag- 
netization of the left barrier is aligned and pinned along 
the z axis, while the other can be freely switched by ap- 
plying a small magnetic field [13] (coercitive fields as low 



as a few tens of Oe can be achieved, e.g., with EuS [181] ) . 
Moreover, we assume the N-layer resistance to be much 
smaller than the FI resistance, so to neglect any spatial 
variation of the spin accumulation within the N region, 
and that tjf is smaller than the spin relaxation length 
in N. Figure 1(b) displays the FI potential profile for 
temperatures (T) above (left) or below (right) the Curie 
temperature (Tcurie)- While in the former case electrons 
with opposite spin experience the same barrier height $0, 
for T < Tcurie the FI conduction band splits leading to 
two barriers with height $ CT = $>o^5E exc , where a =\ ( j) 
for up(down) spin,and 2SE exc is the total splitting of the 
conduction band EH • The spin-dependent transmis- 
sion probability through a FI barrier in WKB approxima- 
tion is given by T a = e~ 2dK " [lj|, where K a = ^2m<& a /h 2 
and m is the free-electron mass. Since each barrier resis- 
tance R<j is spin dependent (R a oc T^ 1 ), RJR] — e 2dSli 
with 8k = Kj, — Kf, so that spin-up current will largely 
exceed the spin-down one depending on SE excl leading to 
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FIG. 1: (Color) (a) Double-barrier tunnel structure consisting 
of two superconductors (S) and a normal metal (N) separated 
by thin ferromagnetic insulators (FI). (b) FI potential profile 
for T > Tcwie (left) and T < T Cur ie (right) (see text), (c) 
Schematic energy-band diagram of the double-barrier struc- 
ture in the parallel (PP) and antiparallel (AP) alignment of 
magnetizations in FIs. A is the superconducting energy gap. 
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a high spin polarization. Yet, in the present double spin 
valve the total current will depend also on the relative di- 
rection of magnetizations in both FIs [see Fig. 1(c)] @. 
In particular, the total current flowing when the magneti- 
zations are aligned in the parallel configuration (PP) will 
largely exceed that flowing in the antiparallel alignment 
(AP). 

In the sequential tunneling regime and within a simple 
relaxation-time approximation we can derive the steady- 
state spin-dependent distribution functions in the N layer 
at finite bias V. In the PP as well as AP configurations 
they are given by 



A^+A/' r + (e M ^r i r s/ )- 1 

Mfl+Afrfr + MT^f)- 1 



f AP (e,V,T) = 



f AP (e,V,T) = 



M+K + iTiTsf)- 1 

e^M+Nr + iTiTsf)- 1 



(1) 
(2) 
(3) 
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where Mi 



Hr) 



= M{e T qV/2,T), Af(e, T) = |Re[(e + 

h)/V( e + *7) 2 ~ A(T) 2 ] I is the normalized BCS density 
of states (DOS) in S, e is the energy measured from the 
condensate chemical potential, A(T) is the temperature- 
dependent superconducting energy gap, 7 is a smearing 
parameter [201 ], and q is the electron charge. Further- 
more, fi (r) = f (e =F qV/2, T), f {e,T) is the Fermi-Dirac 
distribution function, and T 8 f represents the character- 
istic spin relaxation time of quasiparticles in the N re- 
gion. The spin-down tunneling injection rate is defined as 
Y I = (vNRiAtNq 2 )^ 1 , where un is the DOS at the Fermi 



level in N, and A is the structure cross-sectional area [21 1. 
From Eqs. (1-4) it follows that when the injection rate 
greatly exceeds the relaxation rate, i.e., when r^r s / ^> I , 
the distribution functions in the N layer deviate consider- 
ably from fo(e). In particular, in the absence of inelastic 
relaxation (I\r s/ -> 00), f pp (e,V) = f PP (e,V), while 
in the opposite limit (T^T s f — > 0) equilibrium is recov- 
ered. The electronic transport properties of the double 
spin valve are determined by the spin-dependent distri- 
bution functions. In particular, the quasiparticle current 
(e.g., evaluated at the left interface) in both configura- 
tions is given by j pp ( Ap ) — jPP(ap) ^^ere 



ir {Ap \v) 



qR 



deM(e)[fi(e) 



fr (Ap) (e)}, 

(5) 

while the differential conductance is given by 

G PP(AP) { y ) ^ dI PP(AP) /dv _ 

Figure 2(a) shows the differential conductance vs bias 
voltage V in the PP and AP configurations calculated in 
the full nonequilibrium limit (TiT s f — > 00) for d = 4 
nm at T = 0.1T C , where T c = (1.764fc B ) _1 A is the 




FIG. 2: (Color) (a) Nonequilibrium differential conductance 
G vs V in the PP and AP configurations calculated for d — 4 
nm. The inset shows a vs d (see text), (b) Nonequilibrium 
TMR vs V calculated for several barrier widths d. In (a) and 
(b) we set $0 = 0.8 eV, SE exc = 0.2 eV, and T = 0.1T C . 



superconducting critical temperature, Ao is the zero- 
temperature energy gap, and ks is the Boltzmann con- 
stant. In the following we suppose T c <C Tc ur ie, so to ne- 
glect any temperature dependence of SE exc . The curves 
are normalized to the N-state nonequilibrium conduc- 
tance in the AP configuration, G^ p = 2e 2dSK [R l (l + 
e 2dSK )}- 1 , and we set $ = 0.8 eV and 6E exc = 0.2 eV as 
representative parameters for Eu-based FIs In both 
configurations G is almost independent of voltage bias for 
|V| > 2Ao/q (which corresponds to that in the normal 
state), and the ratio G PP /G AP obtains ~ e 2dSK /A for 
e 2dSn ^ l Tnig large difference between G pp and G AP 
stems from the different barrier heights experienced by 
both spin species in the PP and AP configurations [see 
Fig. 1(c)] [9(. By contrast, G is strongly bias-dependent 
for |V| < 2Ao/<7 due to the presence of S electrodes. This 
leads to a tunneling magnetoresistance ratio (TMR) , de- 
fined as TMR = (G pp /G AP ) - 1, which turns out to be 
largely tunable by changing the voltage bias across the 
structure. The nonequilibrium TMR vs V is displayed in 
Fig. 2(b) for T = 0.1T C and for several values of d Q. 
We note first of all that TMR is strongly enhanced as the 
barriers are made thicker, as expected for FI layers . As 
G PP(AP )(0) ~ ( 1 /A )G PP{AP) at low temperature [H, 
TMR(0) ~ TMR at where TMR N (%) ~ 25e 2dSK is the 
N-state TMR. Moreover, while TMR is nearly voltage- 
independent for |V| > 2Ao/q approaching TMRjv, for 
each d it can be tuned and suppressed by several orders 
of magnitude for \V\ < 2A /q. In such a bias voltage 
range, by increasing barrier thickness yields to an en- 
hancement of TMR suppression with respect to TMRjy. 
For instance, with the given parameters, TMR can be 
tuned already over four orders of magnitude for d > 2 
nm, and stems from superconductivity-induced nonequi- 
librium in the N region. The TMR relative variation (a), 
defined asa= [TMR(0) /TMR( A /<?)] - 1, is calculated 
in the inset of Fig. 2(a) vs d, and clearly shows this effect. 
In addition, thicker barriers widen the voltage interval of 
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dent of voltage. Vi ~ 100% can thus be obtained with a 
suitable choice of FIs and barrier widths 



13, 14, 15, 1 



FIG. 3: (Color) (a) Nonequilibrium TMR vs V calculated for 
several values of T. Dash-dotted line is TMRjv. (b) TMR vs 
V calculated for different values of ViT B f at T = 0.1T C . In (a) 
and (b) we set $0 = 0.8 eV, SE exc = 0.2 eV, and d = 4 nm. 



TMR tunability, so that larger d values are to be chosen 
in order to extend the TMR relative modulation with 
V. Note that the finite-bias features appearing in the 
TMR characteristics of Fig. 2(b) reflect how the DOS 
of both S leads contribute to the total conductances at 
finite V through Eqs. (1-4). We stress that nonequilib- 
rium as well as S electrodes are essential elements for the 
observation of this spin valve effect. At equilibrium the 
distribution functions in N would be thermal and spin 
independent, while in the absence of superconductivity 
the nonequilibrium TMR would be virtually independent 
of V (at least for bias voltages of the order of Ao/q). 

The role of temperature is shown in Fig. 3(a) which 
displays the nonequilibrium TMR vs V calculated for 
several values of T and d = 4 nm. TMRat (dash-dotted 
line) is also shown for a comparison. In particular, at the 
lowest temperature (T = 0.01T C ), TMR turns out to be 
tunable with V over about seven orders of magnitude. 
Then, by increasing the temperature leads to a reduc- 
tion of this effect and, for higher T, TMR approaches 
TMRjv which is independent of voltage. The spin valve 
performance is thus enhanced at the lowest temperatures, 
although even at T = 0.8T C the TMR can still be tuned 
over about two orders of magnitude. 

The impact of inelastic relaxation in the N-region is 
displayed in Fig.^b) which shows TMR vs V calculated 
for different rjr g / values at T = 0.1T C . The general effect 
of decreasing r \T a f is to suppress the TMR due to spin 
mixing which destroys the spin imbalance established in 
the N layer. Notably, even in the presence of sizeable 
relaxation, e.g., T]T s f = 10~ 2 , TMR values as large as 
10 4 % and with a tunability over five orders of magnitude 
can, in principle, be achieved. 

The spin-filtering nature of this setup can be quanti- 
fied by the current polarization, defined as Vi = {I pp — 
lf p )/(lf p +l[ p ) = {e 2dSK -l)/{e 2d5K + l). This expres- 
sion holds in both full nonequilibrium and equilibrium 
also in the absence of superconductivity, and is indepen- 



In light of applications, the device could be used for 
the implementation of storage cell elements or logic gates, 
where TMR can be tuned continuously through an ap- 
plied bias voltage. Magnetic field sensors or magnetic- 
field-controlled current switches could be envisioned as 
well. In addition, negligible power dissipation is intrinsic 
to the structure due to the presence of S electrodes, and 
makes this spin valve attractive for low-temperature ap- 
plications. With respect to practical realization, super- 
conducting aluminum (Al) combined with Eu-based FIs 
appear as promising candidates 0, S] , while a conductor 
with a sufficiently long r s f is required for the N region. 
To this end heavily-doped Si, thanks to its extremely 
long spin lifetime 2J, 25 1, seems suitable and could be ex- 
ploited, for instance, in a planar-like structure. We recall 
that r s f as large as ~ 1 /its has been measured in n-doped 
bulk Si (with n ~ 1 x 10 19 cm" 3 ) at 4K fH. We derive 
here a general condition to retain a substantial TMR in 
the double spin valve. By referring to Fig. 3(c), a large 
tunable TMR is achievable provided > 10~ 4 , i.e, 

For the above cited case of Si 
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RiA < 10 A r sf /{u F t N q 2 ). 
with 1/$ - 6 x 10 45 J 



45 T-l m -3 



and 



rSi _ 



10~ 6 s [24( we get 



'W" UAJ - U u 111 <sf 
R^A < {65 ft if) x 10 9 fi/xm 2 for the spin-down specific 
resistance R[A and ijv(nm), a condition that can be met 
with a proper choice of material parameters 

In summary, we have proposed an out-of-equilibrium 
superconducting double spin valve based on ferromag- 
netic insulators. The structure is intrinsically simple 
and can be implemented with present-day technology. 
Besides providing very large and tunable TMR val- 
ues (~ 10 2 . . . 10 6 %) which make it attractive for low- 
temperature magnetodevices, this setup may provide 
physical insight into the spin dynamics at the nanoscale, 
e.g., into the spin lifetime of nonmagnetic conductors. 

We acknowledge F. Taddei for fruitful discussions and 
for critically reading the manuscript, and partial financial 
support from the NanoSciERA "NanoFridge" project. 
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